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ABSTRACT 

The dynamic stability analysis of the airframe for consideration 
in the control problem was made considering six degrees of freedom. 
The intention was to explore a method and determine the extent that 
such a scheme could be used in a system Known to be non-linear at 
larger disturbances. It was determined that a very workable system 
exists whereby differential equations may be simply converted into 
block diagrams and vice versa. The system is workable for various 
flight conditions. It can also handle disturbances up to one tenth 
radian in most cases but is restricted to linear approximations. 

The system is compatible to any multi-loop control problem 
and is not restricted to airframe control. A standardized form now 
exists upon which multi-loop compensation theory research may be 
conducted using either basic hardware or differential equations as 
the basic system. 

The graduate work, for which this thesis is a partial require- 
ment, was performed while the author was attending the U. S. Naval 
Postgraduate School during assignment to the Aeronautical Engineering 


curriculum (Avionics). 
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i Introduction 

The types of motion that result from a disturbance in some 
equilibrium flight condition and the transient motion of the aircraft 
in response to control movements are analyzed in the study of the 
dynamic characteristics of an airplane. Over the course of years the 
effects of dynamic stability in aircraft have been known to designers. 
The mathematical computations were discussed by early aeronautics 
pioneers which included Lanchester, Bryant and Glauert. In spite of the 
fact that dynamic stability effects were known, most people paid little 
concern to this field. As design groups solved other problems of air- 
craft design, the dynamic stability effects were automatically solved as 
a secondary result. However, since World War II, aircraft development 
has pushed into the field of very high performance aircraft. In many 
cases the dynamic stability problem doesn’t solve itself as happened 
previously. Furthermore, the fire control problem has entered the scene 
which has generated a control problem of no small magnitude. 

This work will consider the aircraft in six degrees of freedom and 
explore paths of compensation and solution. It is the prime intention to 
explore the work of Chul) in multi loop servo systems and determine the 
value of it as applied to the aircraft stability. 

The first section will set down a general definition of terms and 
the basic equations of motion as related to aircraft. It will further set 


down the basic assumptions upon which the equations are formed. 
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The equations motion for the airframe which are completely 


derived in Appendix I are: 
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Ww (cos sin @cos WY + sin p sin YW ) (7-13) 
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Consideration of the aerodynamic forces must be resolved into 
the effects they produce on various parts of the body in some flight 
condition and also the resultant moment effects they produce. 

The first or necessary condition is to establish a reference 
condition and proceed from there. A reference condition is some flight 
condition where equilibrium exists in all equations and the summation 
of forces is zero. This is a steady state condition. 

The set of axis upon which this rigid body is oriented is commonly 
called the Eulerian Axis system. General aerodynamics usage modifies 


the z axis so the positive direction points downward. 
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The first rotation moves about the Z axis in a direction WY and is known 
as yaw. The next is the movement about the Y axis whichis term © 
and known as pitch. The third motion about an axis is roll about the xX 
axis and is maneuvered in terms of g . Caution must be used in de- 
termining these positions because the order in using this system is not 
commutative. 

The derivation for the yaw rate R, the roll rate P and the pitch 


rate © are shown in detail in Appendix Two. They are: 


Pp = d -Wsin @ (I-16) 

Q = cos f + W sin g cos © (I-17) 

R = Woos f cos& - S sin Dp (I-18) 
and where small angles are involved 

Pp = DL (I-19) 

oes (1-20) 

Re WY (I-21) 


Similarly, Etkin shows: 


e = Q cos g - R sin L (I-22) 
D =P+Q sin Ptan@+Rceos gF tance (I-23) 
‘ = (Q sin + R cos f )sec @€ (I-24) 


The general equations of motion have to be applied to the aircraft 
motion. The effects of axis movement upon the weight force of the 
aircraft will be shown as well as the reference state condition. It 
must be remembered that the gravity force is always oriented toward 


the center of the earth. The reference body axis may or may not be 
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oriented in this direction. 
From the discussion of the Eulerian angle concept, it is 


determined that a reasonable reference condition would be 


P= Oe xX Aeamensin. Cc ({[~25) 
Xx oO O 

= = i. -26) 
Fy os W cos @ 7 sinQ | (I-26) 
F/=0= Ze cos @, cos Y, (I-27) 


This allows for a steady state condition where the aircraft is not 
accelerating but neither is it necessarily in straight and level flight. 
From the previous illustration of Eulerian angles coupled with the deri- 


vations in Appendix II, the gravity forces can be determined as shown 


below. 
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cosY sin@cos ZG sinYsin@ccs yg cos@cosf "ny, 


= sinY sin fp “10s? sin g 


4 
—_ ; 
_ ~~ 


Discussion and derivation of this general cosine matrix are 
found in several sources (2,3). The overall result of the equations of 


motion is shown below. 
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The above equations are complete except for the external forces 
xX, Y, and Z and the external moments L, M, and N. These forces and 
moments are composed of aerodynamic and propulsion effects and the 
moments due to control surfaces. These equations are very general in 
nature and contain the following assumptions for simplification from the 
most general and complex case. These assumptions are: 

1. The earth is assumed to be fixed in space and the Earth's 

atmosphere is assumed to be fixed with respect to the earth. 


2. The airframe is assumed to be a rigid body. 





3. The mass oi the body is assumed constant during the time 
of the problem. 
4. For an aircrait or missile, the XZ plane is assumed to have 
mirror symmetry. 
The above equations then cover the six degrees of freedom found in an 
aircraft or missile. Coupled with the other auxiliary equations covering, 
roll, pitch and yaw effects along with control motions, solution is pro- 
bably possible. To date, this has been accomplished either by a michine 
using specific coefficients for various terms or several simplifications. 


Discussion of this will be conducted in succeeding sections. 
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23 Consideration of External Forces and Moments on a Body in 

Motion 

In the previous section the equations of motion were discussed. 
The aerodynamic thrust and control effects were simply grouped in X, 
Y, Z, L, Mand N. Itis the aim of this chapter to discuss these in 
detail with respect to the effect on the dynamic control stability problern. 

Consider the term X when the aircraft is in a disturbed condition. 
The reference condition is x and thus X = Xx + AX for some disturbed 
condition. The question is then to determine the A X components. 

The force and moment variables can be expressed in a Taylor's 


S@ties® 


Bio (ar\ur(as\ee-—-- (area a MY 
axfo (gx?) OBI — BE 


For practical reasons the effects of second and higher order terms are 
omitted from further considerations as will be discussed later. Con- 
sidering the motion of the aircraft itself the X forces would be 
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Similarly, 


For convenience all terms will be expressed in letter subscript form. 
It can readily be seen then that considering X, Y, Z, L, M and N and 
the effects of U, V, W, P, Q, and R plus the acceleration terms U, 
ree 

V, W, P, Q, and R there are 72 derivative terms required. This does 
not include the control effects which can add another 36 derivatives. 
Needless to say, considerable effort has been expended in justifying 
simplifying assumptions. In addition to this ,- obtaining reliable co- 
efficients tor various derivative terms is not eee possible. 

It is not the intention of this work to cover the computational 
methods for these various derivatives. The scope is too large. A 
search of all authoritative books 3, 4,5, 6lon the subject yields a 
table of terms which can be reasonably determined. This appears as 
Table I. These considerations are mentioned to bring out the awareness 
of the different causes and effects of external forces. There are con- 
siderable gaps in the table. Determination of these to complete the 
table would be a topic of a complete study in itself. Much work remains 
to be done in this field. 


Furthermore, some of the methods of computing the known 


stability derivatives are questionable. Some computations from theory 
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KNOWN AERODYNAMIC STABILITY DERIVATIVES 


TABLE 2-1 
BlanKs indicate derivatives for which no estimate can be feund. 





prove quite accurate. Others are very poor and empirical methods or 
model testing are used. Even these fail at times and "educated gue: e. 
from previous airplanes are used as a basis of determining stability 
derivatives for a new design. Tail effects are the most notable in thi 
last cateyory. 

In all aerodynamics problems, considerable effort is made to 
non-dimensionalize the equations. The advantages of doing this dre 
of considerable merit. For the most part, itis rather simply done. 

It has been found in both theory and actual practice that 
F=4 Pv*sc 11-4 
F 
Where 5S is a reference area, usually the wing area, V is the absolute 
reference velocity and Cp is some non dimensional coefficient. C.. 
varies with Mach Number and Reynolds Number depending on the force 
_ considered... (lift, drag, etc.) ye is the density of the fluid through 


which the body is passing. Thus a term such as 4w has the following 


relationship: 
Zu' = A of i fares 2 C7 Vs LP V23C, [T= 
ae z : 


Here again V and 5S are the velocity and wing area terms respectively. 
The same concept applies to the moment terms except that 4 moment 
arm term must also appear in the equation. A generalized example is: 


ee TC ft [l=t 
is lo Lin 


In this example M 1s some moment and Cu is a needed dimensionles | 


coefficient to satisfy the equation. jes , V and S are the same as 
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mentioned above while’ is the moment arm. 
Likewise in a stability moment the derivative can e‘‘’* 


form. 
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Another external force to be considered is that of thrust. 
This becomes an important factor in most flight conditions. The 
location of the thrust axis is seldom coincident with any reference 


axis in the body. 
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Assuming the thrust line lies in the xz plane, then the thrust 


forces can be resolved in this generalized manner. 


X =TcosO Il--8 
T T 

Doe =e Sri I-93 
T aj 

M,. =e Ii-10 


For a steady state condition the equations become 


oF = To cos 6 . I-11 
Zon = To sin@ T [i-12 
or = Toe Ti-13 


while in a disturbed condition where the body axis remains fixed in 


the mass, 
= ae 
Xr T, cos @ t li-14 
Zn = Ty sin 6 T IT-15 
Mr = Ty Il=1%6 
and Ty =To+ A [I-17 


Assuming the air density remains constant so the engine thrust will 


not be affected 


IT-18 





These are the considerations for the external forces on a body 
in motion. While the dimensionless coefticient has been used exten~ 
sively by the aeronautical engineering field, it will not be used further 
in this work. The primary mission of this work is to investigaié the 
control response and paths for possible compensation. 


The generalized equations of motion then appear: 
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Thus, it has been shown that the equations of motion have the 
above external forces which will affect the rigid mass in motion. The 
assumptions made are: 

1. The earth is assumed to be fixed in spece and the Eartn's 

atmosphere is assumed to be fixed with respect to the 
Earth. 

Z. The mass is arigid body. 

3. The mass of the body is assumed constant during the time 

of the problem. 

4. For an aircraft or missile, the xz plane is assumed to have 

mirror symmetry. 
These assumptions were made in the first section, In addition the 
following assumptions and simplifications are made in considering the 
aerodynamic forces. 

1. The flow is assumed to be quasi-Steady. 

2. The effects of auxiliary aerodynamic equipment such as 

tlaps, speed brakes, slats and spollers are neglected. 

: 3. The thrust can be affected only by change of speed and/or 
revolutions per minute. No change is considered for chanye 
in fuel mixture, manifold pressure, side slip or propeller 
fin effect. 

These assumptions were made for several reasons. In the first and 
third cases, absence of adequate information on many of the coefficients 


makes ignoring or setting them equal to zero mandatory. In all probanility 
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these coefticients are very small. Assumption two covers special 
cases. While these items may have considerable effect when used, 
they are eliminated from the general case at this time because an itir- 
craft probably uses them less than two percent of its flying time. For 
the special cases, these effects may be quickly added. 

The work in these first two sections is developed primarily for 
background purposes. There is nothing new or startling in this. For 
further examination of this subject, the reader is referred to the works 
in Ref. 7 and 8 which closely parallel the development in this work. 

The aerodynamic force and moment terms carry a primed notation 
for ease in publication. In later sections these terms will be divided 
by mass or moment of inertia. These terms after division will not 
carry the prime (’) notation. There are far more of the latter terms 
used. Thus, it is for ease of publication only and does not affect 


theory in any way. 
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S. Tne Airtrame Equations in a Standard Form 

This section will take the general equations developeu in the 
previous chapters and simplify them. The simplification will extenu 
to an aircraft in straight and level flight and subject to small pertur- 
bations. This is the standard method of dealing with the problem as 
found in present day texts. However, these equations will be used 


in block diagram form and in determinental arrays as developed in 
Chu's work. (1) The results will be compared with present day 
classical methods to determine validity of this determinantal concept. 

The first simplification is the concept of the lateral system 
being independent of the longitudinal system. Therefore, the six 
equations of motion are divided into the two groups. The longi- 
tudinal group includes the directional motion in the X and Z directions 
and rotary motion about Y axis. The lateral group consists of direc 
tional motion along the Y axis and rotary motion about the X and 2 
axis. The axis system referred to here is the body axis systern of 
the airplane. 

Secondly, the aircraft is initially in a steady state condition. 
This is the reference state as mentioned previously. When the cis- 
turbance terms are all zero Equation II-19 reduces to: 
m (Us + Qavc = RoVOi + m= sin © qo = Xo +r T, cos Oa WJ=-1 
However, it was stated that this situation was restricted to lonyi- 
tudinal or lateral systems. Thus for motion along the X body axis, 


a yaw motion Ro and a side motion Vo would be zero. Further it wi:s 
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statec that the aircraft was in a steady state reference condition. 
Therefore us would be zero. Thus equation III-1i reduces to 
m ‘QoWwe) 4 Ma son on = X'o = Ts cos @ T ITI-u 
Sirtilarly, it may be observed that equation II-Z0 for the steady 
state reduces to 
Yo = 0 ill-3 
Likewise equation [1-21 reduces to 
m ! coYo) « 7! +m cos © *To sin@ TeBb6 Tl-4 
The moment equations all reduce to zero provided the rotor effects 
are neglected or considered zero. 
Lo = Mo = No# 0 ITT-5 

Recalling that only small disturbances were to be considered 
and these disturbances were to be in only the symmetrical or lonsi- 
tudinal system for the longitudinal group the derivative due to asy= 
mmetric variables of V, P, R, E and e are zero. Ina similar manner, 
when the asymmetric system effects are considered the symmetric dis- 
turbances are ignored. This leaves the variables V, W, Q, and ¥) 
out of the asymmetric equations. Finally, because of the smali dis- 
turbances, the change angle © will be small so the SDDrOmae On 
cos @ =land sing = PD is valid. 

The result of the previous discussion is subtracting equation: 
Ii-2 from 11-19, eliminating the asymmetric changes and supplying 
the small disturbance considerations. The result is: 
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an = f - Te \ > 2 rv = 
m (nu + Ooa.r Tog + deg) + me @cos@ , = XM + KK + mr 
. : : ‘ i v4 ‘ 
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For the z direction taking ITI-4 from l1-21 and using the same 


assumptions: 
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For the moment in ae y axis, the following applies: 
- oe ag ae fr 3 % ao a 7 ft 
mm = Ball + Page + Mya + Mc} + MAC + va * + Ty o , - 
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To complete the case for the symmetrical group, the action of the 
elevators must be included. From appendix I the elevator hinge 
moment effects are reduced to: 
ee »> es i. ca ‘© fe ‘ Bee 7 é ~ ™ Om 
Te TL 2 Gi + 4 tent Mewd + News @ Meu + Hex? + Ade Ii-3 
: =) Gy] 

In most cases where moment of inertia terms are handled it 
is usually expedient to use an axis system which is the principal 
inertia axis. This means that the product of inertia terms are zero. 
For aerodynamic motion problems, this applies for some cases. () 
This case is not one of them. What is gained in simplification of the 
rotary motion equations by use of principal axis is lost in attempting 
to handle the force equations with the wind velocity considerations 


in Sine and cosine terms. By orienting the X axis into the reference 


wind velocity the result is: 
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where U= ue ra 
Therefore, in the steady state reference condition Vee =U_ and 
LP = we =C,. With the assumption that u, v and w were very small, 


thelr products and squares can be neglected. 
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and since v4 is much greater than 2 Ugu, then VO = Uj =U and thee 
terms are interchangeable in this section. 

As previously stated, the reference situation is a steady state 
motion and includes first order effects only. The orientation of the 
axis oe is @ tinite value while a = Me = Q. Furthermore, the author 
is hard pressed to visualize an aircraft or missile tumbling along with 
a Q, motion although Po exists in spinning missiles and bullets. 


Further x =X = ZS =—~7S XY = = 0 as an assumption. 


g “¢ RPM 


These terms are usually very small. 


The “standard” determinant form as stated in Chapter 3 of 
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Chu’ is: 
+ PP pS Us 7 pe 
| Pro FB o9 is ase! 
ex OE a + Pry Pian ie PE 3° 
oe me te pr 
ao ap —— 32 


where the PF indicates a performance function. The system may heave 
yo 


ne nodes. This example contains only three nodes. The “standard 


block diagram for this form is: 
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Piewre 3x1 STANDARD RLOCK DTAGRAM DRAWN FROM 3x3 DETERMINANT 
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The nodes are labeled in Roman Numerals. The basic rules 


my. 3 


of nodal signals allow several inputs, but only one output. 17 


ns 


hii 
output from the node goes directly to a performance function with a 
single subscript number. Signals may be picked off from the node 
output and fed forward to another node. An example of this is containeu 


in the signal going through PF In this case the signal is also multi- 


oe 
plied by the function PF, 2 and then fed to node Ill. Feedback takes 
place in a similar manner. The signal output from the first system PR, 
is picked off and fed through PPiy back to node I. Similarly the other 
nodes function in the same manner. 

For feedback blocks, the first subscript number denotes the ovi 
put signal source for the feed back and the second subscript number 
the input node to which the feedback signal is going. The notable 
points about the determinant are: 

1. The main diagonal always has terms of one plus the 

direct path performance function times the self-loop 
feedback performance function. 
(1 + PF LPF.) 

2. All terms above and to the right of the main diagonal 

are feedback terms. 

3. All terms below and to the left of the main diagonal ot 

the “standard” determinant form are of a feed forward 


nature. All of these feed forward terms are negutive. 


The term immediately below the main diagonal term is the 














==) = 


> — } aa ' =o -_-_ .«* «=> ©&£ «& 
- — | i » — <<— - a 


— oT”  — —— nv, ez <P 

=> =P > (ee =a 
SS = 

eee 


PP —-— eee 
—_ ieee Om CCUG 








direct path performance function. Thus, where the second 


column main diagonal spot has 1 + PP PP as , the cirect patt 
&: &s &s 


term PF, will appear in the location immediately below \s 


Z 
=—Pr. ° 
L: 
These rules are strictly mechanical. For a rigorous derivation of this 
system the reader is referred directly to Chu'’s work. 
In Chapter 9 of Etkin {4) , there exists a two degree of freedom 
situation where ©and W vary and the speed in the forward direction is 


constant. Tnis is the short period approximation. The equations usin; 


the LaPlace transforms and nondimen ional coefficients are: 
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This author also grouped his acceleration terms such that G, 
na 
included C ‘ and Crap, . Similar as exist for Gay 1 Gig : 


and G . Note that r~ W and ? = yp for this consideratior, Trom 


these equations a block diagram is evolved. 
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Further, Etkin derives: 
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Not it is intended to use Chu's work involving the use of 
“standard” block diagrams and “standard” determinant form to prove 
that the “standard” system applies for the airframe as well 4S u 


multiloop servo system. 
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The Z, M and +j equations of motion with the applied asfump- 


tions ena 2. = 0 are? 
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Recalling that ww = Uo t-n A “and for small angles 
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a) = U, && and O--c, 
The equations in LaPlace notation then become: 
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Ordinarily, it is not necessary to divide by s or s“. How 
ever, it is done here to demonstrate more closely the standurd tor: 

to take the differential equations and arrange them by rows nv 
columns to see the actual physical system is to no avail. The amrunge 
ment does not exist in the above form. It will be shown how to achieve 
this form later. 

Consider the “standard form" array as set forth above. The 


block diagram for this in Chu's "standard form" is: 





Pocnre 383 CUSTLT OF BLOCK DIACRAMING ECUATTONS 
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The above block diagram then is a three node system but the outputs 
from the first and second main path blocks are nothing that is renuily 
recognized in an airframe. Because of the fact that the main bloct: 
following the third node 1s 1, the output hereis © . 

A few simple block diagram manipulations clears the situ ation 
and shows true signal outputs that are physically realized in the :ir- 
frame. Block diagram manipulations allow the changing of a pick ott 
point provided the transfer function is corrected. Further, two block.. 


can be combined into one equivalent transfer function. 
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Adding a dummy node this now changes to: 
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The result of using these two maneuvers rearranges the block a) oj! 


to the follow:ns: 
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RTARRAN ED BLOCK DIATRAM OF FIGTRE 3-3, USTMG MAMITULATIO 
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Maneuvering the above diagram location, the block diagram tukes the 


form of Fig. 3-7, which coincides with Fig. 3-2. 
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It is noted that this block diagram approaches the form set 


down in Eikin except for nondimensionalizing the terms. The dasheu 


area in the diagram is what Etkin labels Gay, and G & y 


‘The determinant from the above dashed portion of the block diagriam 
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Evaluating for the transfer functions © and AX , the results 


oN 
will be in the form GikAiv 
| | | (1) (12). 
This form is the type expressed by Chu and Thaler . the A 
term, which is the whole determinant, reduces to: 
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To consider the 4X term the cofactor 


obtained bv crossing through the first row 


A y of this determinant is 


and the second column. 
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The main feed forward transfer function between the input signal from 


node 2 to the output is +] in this case. 


input to owtput. Hence Gondk 


The transfer function 
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This is also true of node 4 
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It can be noted that this equation III-24 is the same equation as 


IWI-15 given by Etkin save the non-dimensionalizing terms. 
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To evaluate the transfer function & j 7 , the same method is u; ea 
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Here again, excepting the non-dimensionalizing, the result is the 
same. Equation IJI]-26 is equal to equation [II-14 as given by Ltkin. 
The system of "standard" forms then works for the airframe 
differential equations after some manipulation. It can be noted thai 
there is a method of making up this standard form without the necessity 
of the manipulation which will yield physically realizable output sign:s 
When the block diagram of Fig. 3-6 is converted into 4 determinant 
for solution of the transfer function, there appears a Series of alterna-~ 
ting rows where the equation -1+ 1= 0 appears. It also appears that 
these simple equations occur in a sequence where the plus one term is 
placed in the 2n! row and column. nis the number of differenti.il 


equations in the problem. The minus one term is placed in the ?’n! 


rOW w 
and (2n!-i} columns. All feedforward terms (those below and to the 
left of the main diagonal) are moved one column to the right and remain 


in the sime row. This doubles the size of the determinant. However, 


there are several zero terms introduced and the actual work involved in 





evaluation of tte determinant is not markedly increased, 

The next item 1s to look at the lateral motion problem. [* is 
slightly more complex because there are five differential equstions 
of motion instesa of four. 


*) 


Equation IIl-2 with simplifications becomes: 
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The aileron hinge motion equations from Appendix i are 
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And also trom Appendix I, the rudder hinge motion is 
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From these terms a determinantal array can be made. It sneou 


ba 


be notea that this time a set of equations of -1+ 1 = 0 will be insertes: 
in the svstem. This permits direct drawing ot the block diagram and 
omitting the need of manipulation. Table ?-1 then presents this array 


From this Fig. 3-8 is drawn directly. 
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DETERMINANT ARRAY OF LATERAL SYSTEM EQUATIONS FOR SMALL DISTRUBANCES 
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Figure 3-8 BLOCK DIAGRAM OF AIRFRAME 
LATERAL STABILITY SYSTEM AS DRAWN 
DIRECTLY FROM TABLE 3=I 


(Roman Numerals indicate node 
riumber ) 





24e resulting block diagram very clearly parallels thet or the 


U 


A ” . | . : 
lateral system shown by Etkin in Chapter y (4) . Tne “standard” torm te 


achieve this block diagram differs. 
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There is, then, the question of why does this differ from the 
original standard" form set forth by Chu. 

Chapter Iwo of Chul?) Giscusses at some length the idea o! 
arransement of varying performance functions into a “standard” blocr 
diagram. The individual pieces of hardware are at hand with known 
individual performance functions. The orderly arrangement ot these 
permit the construction of a “standard” determinant for the whole system 
From this, a multiloop servo system can be handled in generalized 
results. In the case of the airframe, the determinant formed by the 
differential equations was known and the block diasram had to be 
formed. Admittedly a block diagram (Fig. 3-3) was made up but the 
output signals are not the physically realizable type that are generally 
known to pilots and engineers. The same performance function results 
from the system. Both are mathematically correct. 

The difterence lies in the fact that Chu is taking individual 
pieces of hardware where the functions are known. The blocs are 
made wp and then the “standard” determinant form is made. In the Gace 
of the airplane, the equations of motion are known and s determinant ic 
immediately formed. Each column contains the pertormance functions 
multiplied by the direct function signal output term. Chu's standard 


forms dictate. that all terms above the main diagonal are multiplied by 





the node output or input to the direct function. Steps must be tulen 
then to chuinye this in the differential equations. The addition of tre 
altern:ting -i + | = 0 equations solves this very nicely and in neo wity 
alters tre vilue of the equation. 

A series of evaluations of determinant forms is given in Appeniix 
III. The results show that the addition of these terms in no way chanues 
the final value of a determinant or a cofactor value. 

Keferring to Fig. 3-6, there are other block forms that may chan je 
the values within the blocks, but will not alter the final outcome. Con- 
sider the first node in Fig. 3-6 as an example. The four forms expresse.. 
in Fig. 3-9 are equivalent. The original form is shown in Fin. 3=5(s). 
The self-loop feedback path is combined with the direct path in Fig. 
3-9(b). The 5¢ term is eliminated and the performance functions are then 
in the form shown in Fig. 3-9{c). A further change can be made by 
dividing the nodal input terms by the denominator of the direct path 
feed-forward function. The result is then in the form shown in Fiy. 
3-9(d). The performance functions in each block differ somewhat, but 
the final result remains the same. These manipulations allow the user 
more freedom in the forms of the performance functions. It allows 
chanyes in form to fit the needs of the particular problem. 

The pick-off point for the feed forward terms does differ in 
requirements. It is seen in Figs. 3-9(a) and (d) that the pick-off can 
be either off the output signal of node I or II. The final result will not 


change. However, this is not the case in Figs. 3-9(b) and (c). The 





pict-ofi point must be off the output signal of node II. The tin2! 


result will citter if the output of node I was tapped in the latter c ise 
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Tor converting a set of linesr differential equations into the 

stind .r2 block diagram form, the following rules cn be appliec. 

1. Between each differential equation, an equation =! + 4 =u 
shall be inserted. An additional equation -1 + 1 = 0 snail 
be added at the bottom of the set so the number of addi 
tional equations shall equal the number of differential 


equations originally in the problem. 


> 


The arrangement of the values of these added equations in 
the determinant shall be such that the positive one value 
is always in the main diagonal block. The minus one term 
shall always be placed in the block immediately to the lett 
of the main diagonal term. 


Ali other blocks in the rows of these added equations shill 


Cad 


be of zero value. 

4. Ali finite performance functions of the original differenti:! 
equations which are located below the main diagonal shali 
be moved one place to the right. The vacated blocks shali 
be replaced with a zero value. This last rule is options! 
to the user. It does save the possibility of error whicr 
could appear because of cases shown in Fig. 3-9(b) and 
(d) . 

With the simplified rules of making a standard block diagrim 


from the standard determinant, a valuable tool is developed for muny 


purposes. I[t cain provide a student with a means to draw and visualize 
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a compiex feedback system such as an airframe from a group of differen 
tial equitions. For the research groups, it makes analog computer 
programming a simplified matter with respect to amplifier uSage. in 
this latter crise, pitfalls do exist. The airframe is no simple mitter 
The reacer is referred to Howe is) for an excellent coverage of the ‘ir 
frame problem for analog computers. Most important, though, is the 
fact that this standard form of Chu’s is valid for any multi-loop linear 
feedback system whether it is electrical, mechanical or a combination 


of the two. This then puts a regular form to the multi-loop contro] 


System trom which design or synthesis can proceed in an orderly manner 


4zZ 








4. The Control Problem for an Airframe 

The previous three sections have discussed only the ‘1irfr «me 
itself. It is quite reasonable to say that while these point out miny 
of the responses, they achieve little toward solving the control proble~ 
In simplified terms, the control problem is resolved into the block itu 


gram of Fig. 4-1. 
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Indeed this diagram is an oversimplification of the aircraft 
control problem. However, it does point up the fact that the previous 
three chapters have dealt only with the controlled system block. Any 
complete study of the airframe control problem must include the other 
block components of the amplifier and the error detector. The probier. 
is a multiloop type and of no small magnitude. 

For the control problem, the diagramming becomes a simple 
extension of the block diagram developed in Section III. Consider the 


Standard determinant for the longitudinal system. 
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follows 
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The determinantal array for the longitudinal system is » 
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From this a block diagram can be made up from the determineni i} 
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This diagram i's of the three degrees of freedom longitudinal 
case. Combining this to a control system, the following diagram 


results. 
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Thus. the control problem for the simple attitude O is 
presente in Fig. 4-2 and 3. The two represent the different point o} 
view of the aeronautical and the electrical engineer although the wy ter 
is the same. In Fig. 4-2 the aeronautical engineer is faced with the 
primary problem of compensation within the airframe. For the electrical 
engineer who inherits the fixed airframe, the attention would shitt to 
Fig. 4-3. In this figure, Gy would represent the pilot response in mo ing 
the control lever. Gy would represent the servo system response of « 
hydraulic or electric system. G2 would represent the error detecticn 
and correction signal gained from an electronic element of some Fina 
Similar situations of block diagramming would exist for the change. in 
speed and altitude. For the most part, the compensation schemes to 
achieve stability or improve stability are fixed with respect to paths. 
The change in aerodynamic coefficients will vary the paths in the air- 
frame response group. If the electronic devices are to be used, then 
the pick offs must be at the system outputs and the signuls fed bac 
to near or at the input or reference signal node point. 

In Fig. 4-3 the pilot constitutes the outermost loop. However, 
this is not mandatory. Ina fully automatic system, the pilot circuit 
is opened. Few paths exist where both aerodynamic or electronic 
means sate be used to achieve stable control. One can be noted on 
Fig. 4-2. That is elevator position response 7 which could be altered 
aerodynamically by changing parts of the present feedback block or 


adding another electronic control block. It 1s then fed into the force 
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input noae marked ZA Fe. There is no guarantee thaz this path solve. 
any stability problem albeit the flutter problem may be Aided. I[t merely 
serves 2S an example of two possible types of compensation over the 
same path. 

The lateral systems are not shown or discussed in this chapter 
but the same principles of control systems apply. It should be remem - 
bered that the airframe response shown here covers the small pertur- 
bation response concept and thus is greatly simplified. The following 
sections will deal with expansion of the concept for some workable 


systems. 
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a4 Lxpansion of the Standard Form for Greater Disturbances 

The system developed by Chu was for lineir functions. It 13 
the purpose of this section to explore the equations of motion ing 
determine whether the equations of motion can be used for linger cis 
turbances than considered in Section III]. If so, then to what dqesree 
can this system be used. 

Tor all practical purposes, the small angle approximation held: 
quite well up to one-tenth radian. Therefore,& ~ sin © 
cos@=~1. With this simplification in mind, Equation [I-19 is then 


rewritten. 
m (1 + OW =o RV) + me sin ©, = Pe (cos@y sin.) + 
g 
2 mr (cos, cos aXe SG eee + 


Ei, reer +7. — ie cos U+tT Pes 
4 OT * * cpm & rey C8 er 


The thrust forces are combined into the forward velocity motion 
u and a steady state reference condition of 


i bb m RV *me si = X +fT7 sco Vat 
rn ( on ey 2 sane. 2 5 °° 7 — 


is assumed. 


Subtracting V-Z from V-1, the following equation results 
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The product term qw is dropped as a small value procuct situ ition 
The same epplies to the vr term. 
Converting to La Place notation and dividing by the mass the 


followin3 eyuation will result 
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The y and z direction force equations will form up in the same manner 
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The moment equations can be formed in much the same manner From 
Section Two, Equation II-22 would convert to this equation: 


op AR we i SP... 5. 
B f7-P) (sg of cP? @) = i (sP, @ «= 5 o¥)-— € ¥ 


/ aoe fi 
} 7 


as Wi 
- (sl * Le G - Ing ~ ‘shp) - (sh € oa 7% 


In a similar manner Equations II-23 and II-24 can be brought to a 


similar maximum linear condition. They are shown in the following 


equations: 
- 8 
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Also the control hinge moment equations III-9, III-31 and 


III-32 are recalled in their same form. 
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These nine equations represent the linear equations of motion 
of the airframe and include all cross~coupling effects. These eqixation. 
still retain the requirement of a linear approximation. The equation; 
with the author's modification of adding the equations -l+1]1=0 mute 
up Table ¥V-I. This is an 18 x 18 determinant. From this table the 


block ciagram of Fig. 5-1 can be drawn. 
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TABLE S-I DETERMINANT ARRAY OF AIRFRAME WITH LINEAR APPROXIMATIONS For six Decrees °F FREEDQ{. = a a . ole i . 
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FIGURE 5-1 BLOCK DIAGRAM OF THE AIRFRAME AND CONTROL SYSTEM USING 


THE STANDARD FORM 
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2hts plock diagram then represents the complete ttrir ue 
response muintiining linear approximations. Lven then this ysiter 
has limitations. As an example, consider some airplane tiying with 
a steady state wing angle of attack of 15 degrees. (assume 2 “con- 
vention” wing and the reference axis is the zero lift line) Any pilot 
will quickly testify that an increase of 5 degrees angle of attack. at 
that point will probably result in a rather wildly oscillating mineuver 
popularly known as a spin. This condition is generally outsice of the 
range of the aerodynamic forces and moments usually cormputed. Some 
aerodynamic coefficients such as lift (Cz) are decidedly non~-line=r in 
this region. Therefore, the restrictions imposed on this derivation of 
terms must include a restriction that the airframe is not operating «+1 


the very edge of its flight envelope. 


& 


The system has been diagrammed directly from Table V-]. Tr 
terms in Column One are all connected at node I of Fig. 5-1. Likewise, 
the terms in Column Two of Table V-i are all connected to nowe Tl. cxtra 
nodes appear in the beginning and end of the block diagram. These ure 
inserted as provisions for control functions. These control nodes are 
lettered alphabetically. The block diagram contains the required loop 
for aileron and elevator control. The rudder control loop can be esily 
inserted when desired. 

LXamination of Fig. 5-1 reveals many interesting points here~ 
tofore knuwn but sometimes only by intuitive reasoning. The concept 


of separate systems for the symmetric axis and asymmetric axis 





yster7s is very well shown. The level fliyht steady stite condition 


uSeu in Section ihree stated trat 


cxarining the paths between nodes VIII and [¥ reveal tnt all 
of the performance functions are zero under these specilied condition. 
of steacy state flight. There is no direct, feedforware or teecdbac: 
paths between nodes VIII and IX. Thus the restricted case of stenudyv 
state level flight conditions with small perturbation is clearly divided 
into tie two systems The longitudinal system includes the forwir: 
velocity, vertical velocity, pitching motion and elevator position. 

he block diagram for the system can be further restricted to two 

Gegrees of freedom by holding the X velocity change, u, at zero. 
This allows change in vertical motion and pitch. This would then be 
Figs. $-6 and 3-7 of Section III. 

The so-called lateral airframe response system comprises the 
nodes trom LX to XVII. The rudder and aileron motion, siceslip, y iw 


and roll comprise the motion for this sytem. With the steady state 


Wi 


level flight conditions previously mentioned, the blocr. diagram i 
identical to Fig. 3-8 of Section ITI. 

An interesting note of these separate systems is the crose 
coupling terms thet begin to arise with larger disturbances. An ir 
portant one is the pitch due to sideslip. The effects sire labeles me 
in the diagram. Several aircraft are subject to this. Usually, this 


is 2 tucking or nose down motion that appears as cideclip angle ( 4 


is 
Con 





increaces. In the simplified systems this efiect was neglected 11 
is of nealiganle value for very smell disturbances, but Cun Leco”e 
noticezhdie by the time the sideslip angle Cs ) reaches five leuyree 
vne effect exists and is shown by this expanded system. 

=uxcellent methods of analysis are available for these standurd 
forms brought forth in this work. Methods of analyzing outputs in 
response to one or more inputs are available. Work by Thaler! con: 


tains < very useable system. 


Where iv ze the output signal at the n th node 
A 4 characteristic determinant 
do i the cotactor of the determinant with °s" 
n the node location receiving the input 
and "n" the node location of the output 
signal. 
& | | 
Gi, = the performance function on the input signa] 


prior to entering the input node. 


tof 
\ 6B 


the input signal 


i > 


the direct path performance function between 
the output node and the output signal. 


Thus 4 single input signal would yield an output signal 


; + a’ 
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For two or more input signals additional subscripts would be used ind 
the tunctions are additive. The rules of superposition of linear outputs 


justify this method. 
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TRIS @uatlon would give the resulting output of one or More inputs 


it various nodés in the system. 
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The Jeygrees of freedom can be restricted also. In the bloc}. 


diagram, (Pig. 6-2) consider a cut made in the path at point “‘s“. 
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it 5 Opening in the path now eliminates the sign] to tne utrect 
pits, the feedbick and the feedforward path. The result of tris "n*ve 
determinant is the insertion of zeros in all] the performance functions 
in ice colurn. It leaves only the constant, one, in the main dtagon :! 


block. The reader will recall that the main diagonal has the term 


eo this is shown below as opening the output path from , 


a 
one. 
<3 0) 


lLuavPeAt 


third node. 


x x 0 x x x 
The output path from node three is opened at point s and the ceter- 
minant is changed in this manner. 
A special case occurs for the diagram illustrated in Fig. 
The second node has only one input signal. Because the path was cut 


a0 


at point "s", the output from the second node will aiso be zero 7 


mex S 


all feedforward values from this node will be zero. The cirect anc 


it 


a 


feedback paths of the second node are non-existant to begin with ind 


remain so. Tre determinant form would appear as now shcwn. 


Ct 
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x x 0 0 Xe mx 

xe. x 0 0 x ex 
Simple manipulation now reduces the six by six determinant to a 
four by four determinant. Therefore, when it is desired to consider 
a problem with one or more degrees of freedom held at zero, the simpie 
method of opening the output signal path from the node and placing a 
zero for all the performance functions in the column leaving only the 
constant value one in the main diagonal point accomplishes this. 

Referring to Section II, the reader can observe the general 

equations of motion for mass particles and examine one such as 


Equation II-22 which is repeated below. 
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Note the third term QR is a produce of (Q, + gq) (R «yr ). The pro- 
O 

duct when multiplied is QR, 4 Qer + Ro + qr. From the equations 

of GBesmien I: 
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When the srnall angle approximation is no longer valid, the Lhandlin: 
of this product term becomes rather involved. Furthermore, the ero 
Gynamic coefficients of the Lee Ly. L etc terms lose their linexr 
approxirations over large changes in flight conditions. Most aero- 
dynamic coefficients are of a decidedly non-linear nature over the 

full span of an aircraft's maneuvering capability. The coefficients 

are definitely non-linear with Mach number when the velocity ringe 
covers the subsonic, transonic and supersonic spectrum. The reacer 
can observe the general equations of motion in Section II. Considering 
the complexities of incorporating equations V-14 and V-15 for large 
movements «and the non linearities of the aerodynamic coefficients . 

it is not long before one realizes that the non linearities are numerous 
and cover all the equations. The prospect of block diagramming for 
the purpose of placing all non linearities in one block to be handled by 
a describiny function scheme is hopeless. 


The possibility of doing the problem with an analog comp.u.ter 


and using no approximations is fair. Undoubtedly the task, if attempted 


will be of considerable magnitude. 

Thus, itis seen that this system of standard determinant arrays 
and standard block diagrams can be used to handle the airframe probler 
as long as a linear approximation is maintained. The system is flexibie 


and can cirry varied conditions of multiple inputs. It can also be notes 
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that for the six degrees of freedom problem over a ]arge variation, 


it is of limited value except for diagramming a computer problem . 
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G | Prospect and Ketrospect 

The airframe motion problem has been derived quite zene: -llv 
in the first two sections of this work. The attempt was made to maintain 
as many effects as possible in the general equations of motion. ‘%:is 
was Gone so thut the basic equations of motion could be used in « 
variety of airframes varying from missiles to ground effect machines. 
In the succeedirg sections the author has confined his discussions to 
equations of basic airplane types. Rotor effects were ieee and 
thrust effects were combined with drag to determine the effects along 
the X axis system. This should not imply that the equations of motion 
cannot be used ior other airframe types. They can in all cases where 
the assumptions of a rigid body, constant mass, etc. are still valid. 

ihe author has several thoughts on areas of further exploration. 
He will discuss them briefly and mention any work of possible interest 
in connection with this field. 

The primary aim for control engineering in this field is the 
development of an orderly process for multiloop compensation. It ics 
shown clearly in Fig. 5-1 that several paths of possible compensStion 
exist. What path or paths are the most successful and/or easiest to 
achieve multiloop compensation? Some work has been done in this 
field. Specifically, an electronic compensation of a missile using thi, 


determinant method was done by Anderson and Roane. il 


They contines 
their etforts to the roll and yaw coupling problem for a specific vehicie. 


Their efforts achieved a satisfactory solution to meeting the missile 


C3 





specifications, but no definite method of analysis as to the hest or 
poorest prth to compensate 2 multiloop system was generstei. A 
possible approzch would be an examination of aerodynamic coefficients 
and their known effects. In chapter ten of Perkins and Hiee? there is 
a discussion of the fact that drag or coefficient of drag markedly iaffecr.. 
the damping in the longitudinal phugoid motion of an airplane. The 
airfrare short period is affected largely by the elevator hinge moment 
coefficients C and © . An examination of these factors anc 
DA oF a 

others might lead to a pattern which could successfully soive the 
problem for any multiloop system. 

The block diagram development should provide a basis for 
analog simulation of the linear problem expressed in sections III ana 
V. The analog computer arrangement may be further expanded to inclcie 
some of the non-linear effects using function generators. This proplen. . 
allowing six degrees of freedom, would overtax the capacity of the 
analog equipment presently available at this institution. If an analog 
computer with a 100 amplifier capacity and 50 function generators was 
provided, the problem could be nicely handled. Most assuredly, com- 
binations can be made by grouping various functions which will cut the 
above numbers by one half. However, various outputs will not be 
available for analysis and many signals which are not physically 
realizable in the actual system will appear as outputs. 

The high speed digital computers could be utilized in a project 


of statistical evaluation of compensation schemes. The computer's 











_-_ — —  — om - 

TE — => & 

es > = = gi > 
Ee += > 


i = 46. “> — 
OO — 
oe emma 8 Ce em 
> eee oe ae 
ow = 4b a «s-< 
Se —_— 
ee ee ee 
eee ee a tt & Ol 
Set eet ie 
—i- << = 


» 








i 





cap.tcility of producing hundreus of solutions per minute woulu aliow 
ex.r.in ition of cata to determine trends and patterns for vurivu. co” 
pens::tor schemes. A successful compensation criteria might ue 
derive} in this manner. 

AN serodynamics group may examine four items of interest in 
this fiel?. The control surface hinge moment equations used in tnis 
work were extremely simplified. Exploration of these effects to 
attain more refinement woul?! be justified. 

in section II, the author pointed out the absence of compu 
tational methods for various stability derivatives. FPurther, the 
assuription of a quasi-steady state condition is made in corputing 
the known cerivatives. A study by etkin (2) cites an improved syste 
to attsin better accuracy of the aerodynamic stability derivates. 
Further work in this area appears justified. 

There is also a need to determine the effects of aeroelusticit, 
on the problem. To this end, instrumentation and flight testing ‘in 
airframe might prove useful. The capability of the Naval Air Facility 
to maintain and operate many types of military aircraft should permit 
use of some airplane or aircraft for this work. An additional use woiila 
checking the effects of automatic stabilization equipment in use on the 
airframe, The HSS-1 helicopter has an excellent control system for 
this type of research. However, the equations of motion used in this 
work will require some refinement to control the varying inertia problem 


caused by the tilting rotor disk. 





A final topic for consideration is that of signal flow theory 
This ides, which has been advanced in recent years, may be com 
patible to the same “standardized” system as presented here for tre 


block diagram method. 
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Conclusion 

it has oeen conclusively shown that Chu's “standuird” bloct 
diagram and ceterminant method can be applied to a linear system 
other than @ servo control system. With some modification, any séi 
of linear differential equations can be block diagrammed into 4 stand ra 
form for general analysis. A very workable system exists for analysis 
of a multi-loop control system using one or more input signals, The 
airframe stability problem can be examined for the linear approximation 


using this very effective system. 
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AFFENDIZ 1 
BASIC EQUATIONS OF MOTION 

sve busie concept. of mass in motion and the equxtions ci 
motion cain be derived in vector notation. This is done usin3y mas. 
particles in motion and deterrininy their effects in two separate axis 
systema. The basic reference system shown in Fig. Al-l is the *turre 
axis. This is the non-rotating reference system and is also known is 
the Galllean Axis. The axis system is fixed in the rigid mass but may 
rotate about the fixed space axis system. Itis the relative system. 
A special set of notation is used here in this appendix which differs 
from the body of this work. 


the dcetinition of the special terms in this «appendix is as 


follows. 
Wg A, . . 
i = the vectorial unit along the X axis 
—- &, | 
j = the vectorial unit along the Y axis 
<o A, 
k = the vectorial unit along the Z axis 
“Fy, A, ¢ . 
i = the vectorial unit along the X axis 
> A se 
j = the vectorial unit along the Y' axis 
kK = the vectorial unit along the Z axis 
& ine . 
0 = the origin of the relative axis system 
* & - . 
0 = the origin of the space axis system 
*e bs | | ro 
R = the vectorial distance from the space axis oriyjin to the 
relative axis origin 
—> LQ | 
i == the vectoria! distance from the space axis origin to 4 


mass particle in space 
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I] Be 


the vectorial distance from the relative axis origin to 
the mass particle 


“C4 


A 
dm = the mass particle 
ie A 
P. = the vectorial distance from the mass center to the origin 
oi the relative axis system 
—_—_ AN 
{) = rotational velocity vector of the relative axis system about 
the space axis; an absolute velocity vector 
—? A 
WwW - absolute rotational velocity vector of a rotor whose axis is 
attached to XYZ axis system 
> & 
ie, = momentum vector for the entire mass body 
> A . 
H = moment of momentum vector with respect to the space axis 
> A | 
G = moment vector with respect to the relative axis 
—> A 
iF = the summation of the external forces on a body 
4 


ft the summation of the internal forces on a body 

This appendix will derive the basic equations of motion for .iny 
mass moving in space. After the basic equations have been developed, 
the effects of control actions will be added. The final equations wil! 
then be converted to standard aeronautical terminology which is used in 
the body of the work for aircraft stability considerations. 

These equations have no doubt been derived several times in the 
past one hundred years. However, the author was unable to find any 
clear derivations in any of the reference material covered. Therefore, 


the derivations in this Appendix are submitted for the reader's interest 


and for background information. 
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For an airframe with a large rotor (ie helicopter, STOL aircraft or GEM) 
consider addition of the rotor to the airframe using the mass center of 


the aircraft oe the cia of the axis plage 
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/ Fioure AI-2 VECTOR DIAGRAM OF 
i ROTOR MOTION 


P« is fixed in the rigid body, The rotor spins at an absoluts 


angular velocity of but the mass center of the rotor system is at 
‘> 
~~ 


my 


point P, Therefore, Ff andM, = 0. It is also assumed that 
— 
OO, = constant. Therefore, the additional terma.x(ax/,) 
enters the problem, 
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——j> 5 
[@ Fe fe 
~~: ia Z 


is a basit identity, therefore 
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Note: coordinates i, j, k are those in the body and not the inertial 
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or earth system. The values of—f.are primed for the same 


reason. It is quite pessible that fo will be 0 considerin:; 
y: 


alignment of the rotor on the xz plane. 
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In resolving this, it must be noted that the mass centers differ 
from those of the combined masses. 

It should be noted that the equations for the acceleration of 4 
particle also contain rotational and relative velocity terms, however, 
P and Pp = 0 as these particles are in a rigid body. Further, the 
summation of particles about the mass center makes 
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Moment terms may be considered using Moment of Momentum. ; 
From this the differential may be taken and the moment is then 

derived. Consider H to be the moment of momentum of the system with 


respect to the inertial (non~rotating axis). The derivation will include 


a rotor spinning on the rigid body. 
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Figure AI-3 MASS DIAGRA.. 
FOR DOMUNT ANALYSIS 
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The angular motion of a rotor in motion, but attached to the body will be 
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Henceforth this relative term will be noted as K (KAPPA). 
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The term under the integral will be known as H 
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A further assumption is stated at this point. This is that the 


point of rotation of the rotor is fixed on the body and remains so. The 


shaft does not have nutation about some axis but remains fixed on 3 
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mount of the rigid body. Therefore, the le of the rotating body is 


zero. Also, logically, in a rigid body the mass does not shift hence 


ye is also zero. 
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In crossing a vector with itself Vo x Vo the result is zero. 
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Furthermore, crossing two dissimilar vectors in opposite order yields 
opposite signs and equal magnitude. 
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Incorporating these two facts, the equation for G can be simplified to 
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For any ordinary aircraft or missile, the product of inertia terms Ixy 
and Iyz can be considered zero because of the mirror symmetry aspect. 
Ixz may be zero if the principal axis is selected. Thus with no loss in 


generality, consider Ixy = Ixz = Iyz = 0 
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However, including the effects of Ixz, the product of inertia term, 


Using standard aerodynamic terminology, the following items are 
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These equations concur with Duncan, Page 66, and Etkin, Chapter 4 
rs 5 an e =p 
1o consider the rotor effects it should be remembered that ww. 
LL+j< where K is the relative motion of a rotor on the body. Thus 
for a rotor there will be a combination of terms considering the absolute 
angular velocity of the OAYZ system plus the angular velocity of the 


rotor with respect to the body. 
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It should be noted, that for a symmetrical rotor or disk the 
product of inertia terms are zero where the disk is at least a three 
bladed propeller. The two bladed propeller has a product of inertia 
dependent on its angular position. This case will not be considered, 
but the three or more bladed propeller situation will be used in all 
considerations. Hence with a dynamically and statically balanced 
rotor fe = 0 and the product of inertia is zero. 

The moment of inertia of the disk must be considered from tte 
axis of the rigid body. Therefore, the moment of inertia must include 
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The moment of inertia terms for the rotors will be thus defined in the 


most general sense 


A . a : 
Ix rotor = Ix + md ‘ and aligned with the reference axis 
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b = ly rotor 8 Iy + md*, and aligned with the reference axis 
1. a a | i. 

c = iz rotor = 1z+ md, and aligned with the reference axis 
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Recalling the derivation of moments for the rigid body without the 


rotor, the same method applies. 
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Combining this to a system consisting of a rigid body plus a rotor, 


the following results occur 
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Combining and using aerodynamic notation. 
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This result concurs with Etkin, page 116. 

Therefore, using a standard aerodynamic notation the moment: 
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and the scalar quantities are as follows 
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thus, the force and moment equations derived in this appendix 
are applicable to a rigid body of constant mass where the hody is in 
motion. it also considers the eftects of a three bladed or more rotor 
system which is attached on some fixed point ot the rigid body. The 
rotor etfects may be included cr discarded depending on the magnitude 
of rotor inertia and the relative velocity of the rotor. 

it should be noted that the derivations are done in @ general 
form, but in the concluding eguations, the usual aerodynamic notation 
is substituted to bring the equations to this specific problem. Further- 
more, miscile or aircraft shape and density has been considered sym- 
metrical about the xz plane. Thus, the product of inertia terms Ixy and 
Iyz have been considered zero. No generality is lost by this and if, 
for some reason, the body is sufficiently unsymmetrical, the terms cin 
be quickly added in. 
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The LaGrange equation of motion in a moving reterence trame is 
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APPENDIX {1 
ROLL, YAW AND PITCH 

The derivation of the yaw, roll and pitch rates is given here to 
show the actual motion taking place. Frequently the motion is con- 
sidered small and small angle approximations are used. These will be 
noted at the conclusion. 

First, consider the Eulerian Axis used in the aeronautical field 
shown below. The rotation takes place in the order of yaw ( ¥ ), 


pitch ( © ), and roll { D ), 





Figers AII-1  LCULGRIAN AXIS SYSTEM SHOWING RCTATION ABOUT THE 
Z, Y, AND X AXIS 
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The inttial conditions 1 Ky: * ana Zo and the relations of 
—~. ede anu DL c3in be related to P, Q and R of the body axis at 
xX, Y, Z inthe above diagram consider the angular velocity P whict 
is about the OX axis. It must be the sum of the velocities about the 
OZ): O”* , ana OX. axis respectively. 


2 , TH | 
OZ, makes an angle ~> + @ with OX. The rate thus 


W ens (= a G)e- Wein G. ma is perpendicular with OX so 


y Le ; : 

© cos it = 0 and OX is parallel with OX or the motion 
eee 2 

Oo ot Whe SF 


Thus the summation P= & - Y sin @- 


To evaluate Q the angle YO4Z, must be determined. Specifically the 
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yy 2v | ; . 
value ¥ cos 104. is desired. 
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Thus the contribution of O04, axis motion is yp sin é COs Gs 


The angle between OY and OY is p so the motion here results as 
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© cos gL . OY is perpendicular to OX, hence the angular motion 


7: cos iL = (0, 
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Thus Q = & cos p “ p sin D cos © 

In the third angulir velocity motion R the yaw motion appears as 
y cos cos © while the motion of pitch is 6 Cos (x r g)=-6 sn Fe 
The rolling term vanishes D cos = = Q. 
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Thus R=- © sin g + W cos f cos © 





Hence the relationship of the body yaw, pitch and roll rates 


(P, Q, «nd RK) to the designated reference system angular velocities 
4 , & and D is: 
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p= O- P sin @ 
Q = Ocos H+ Y sin @ cos © 
R = cos G cos@ - G sin Y 


These results concur with Duncan page 81 and Etkin page ilt¢ 
it should be noted that for small angles the approximation below 


is often usec: 


P= ¢ 
qr e 
R= ¥ 


A second consideration which must be made regarding movement 
about the reference axis is the velocity or displacement with respect to 
time. It is desired to obtain the velocity with respect to the fixed frare 

f 


or inertial axis i* x, j*y and ~*2. 


Consider first that in scalar values 
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It should be noted, however, that in the three position system shown 


on the previous paye that 
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and U. =U,cos W -V,, sin ¥ 


Combining the |, 2 and 3 position terms 
dk® =U cos. © cos W +7 (sing sin@ cos W - cos D sin ¥ ) 
+ W (cos sin ©@ cos W + sin@ sin ¥ ) 
dy* =U cos © ain Y+V (sin D sine sin Y tr cos f cos W ) 


+ W (cosZ sin@ sin W - sinY cos ¥ ) 





dz* = -U sin@ + V sing cos@ +W cos D cos © 


Thus, the velocity or rate of change with respect to the inertial 
or earth's reference axis in this cuse is expressed in these terms as 
derived. These velocity terms sre derived in Etkin Page 102 and coincide 
with the derivation shown. 

Therefore, this appendix accounts for the changes in angular and 
directional displicement of 2 rigid body with respect to 2 fixed reference 


axis system. 
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APPENDIX 11] 
NOTES ON DET: RMINANT MANIPULATION 
The bocy of this project contiins several determinant expinsions 
where the orcer is doubled. These manipulations do not alter the vilues 
of the determinant and the resulting values of the performance functions 
are the same. 


The simple determiniunt 


Al fl Cl 
Az 33 ol= A 
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evaluated will yield the following result 
ALBZC? + BICZAS + CIAZBS 
A A3=! 
= {> CiB2ZA2 - AiC2B3 - BIAZC3 
This expansion by the ecuations -1 + ] = 0 will result in the following 


determinant. 
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The resulting evaluation will yield the same value as result {A3-1) 


Furthery ore, if the terms below the main diagonal are allowed to 
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risht as is .one in the main body, the resulting A is still the 3 ure 


ane results of evuluating cofactors tor particular performance 
in uncnan ec. 


functions rem: For eximple it is desired to evaluate the 


cofactor “é\ _. of the original determinant. Thus: 
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The result of Equation A3-3 equals Equation A38-2. With the values below 


the main diagonal! left in their original columns, the following occurs: 
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By inspection of columns two and six, two reductions are made to the 


following: 
=(-1) |aA? 6 C2 
A. (-1) JA 
0 i 0 
| AS t ex 
A 2 CU /A2CS - A3CB / A3-4 


The resulting equation (A3-4) is the same as equation (A3-2). 

A similar examination using a four by four determinant system 
was mace. This was done to insure that the simplified effects of a three 
by three or lower order system did not cover up an error in = four by four 
or higher oraer determinant. The same manipulations can be performed 
and the resultant value of the determinant is not altered. 

Thus, the individual has the option of placing his values in 
either column. The location of the values is dependent on desired 
location of the pick-off point on the block diagram. The final result 


is not altered. 
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